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Quantum weak energy inequalities (QWEI) provide state-independent lower bounds on averages 
of the renormalised energy density of a quantum field. We derive QWEIs for the electromagnetic 
and massive spin-one fields in globally hyperbolic spacetimes whose Cauchy surfaces are compact 
and have trivial first homology group. These inequalities provide lower bounds on weighted averages 
of the renormalized energy density as "measured" along an arbitrary timelike trajectory, and are 
valid for arbitrary Hadamard states of the spin-one fields. The QWEI bound takes a particularly 
simple form for averaging along static trajectories in ultrastatic spacetimes; as specific examples we 
consider Minkowski space [in which case the topological restrictions may be dispensed with] and the 
static Einstein universe. 

ff") , A significant part of the paper is devoted to the definition and properties of Hadamard states of 

spin-one fields in curved spacetimes, particularly with regard to their microlocal behaviour. 
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I. INTRODUCTION 



In common with all observed forms of classical matter, the spin-one fields described by the Maxwell and Proca 
equations obey the weak energy condition (WEC). That is, the stress-energy tensor T a b obeys T a i,v a v b > for all 
timelike vectors v a . In classical general relativity, energy conditions such as the WEC play a key role in many 
important results, notably the singularity theorems of Penrose and Hawking Q. Moreover, since any metric 
t— ( . solves the Einstein equations for some choice of stress-energy tensor, it is arguable that general relativity has limited 
C*~) ' predictive power in the absence of such conditions. 

However, all classical energy conditions are violated by quantum fields, as has been known for many years 
£^ ■ Typically, the energy density at a given spacetime point is unbounded from below as a function of the state. Specific 
examples of negative energy states are provided by highly squeezed states of light in quantum optics, the Casimir 
vacuum state for a quantized field between uncharged perfectly conducting parallel planar plates and the Boulware 
vacuum state outside a black holes. Replacing the classical stress-energy tensor on the right hand side of Einstein's 
5h i equation with its quantum expectation value, we must therefore allow for negative energy sources, raising the pos- 
sibility of exotic phenomena within the realm of semiclassical gravity. For example, negative energy can be used to 
J> , maintain static traversable wormholes 0, § , create naked singularities @, Q , travel faster than light 0, || , and travel 
k>( ' backward through time 0, fTojl . In negative energy models studied by Parker and Fulling it is even possible to avoid 
the cosmological singularity [lj, thus threatening to overthrow the classical singularity theorems. One might also 
expect that negative energy fluxes could lead to macroscopic violations of the second law of thermodynamics ^2 • 

Of course, macroscopic violations of the second law are not observed in nature. Motivated by these thermodynamic 
considerations, Ford [Tj] deduced that the negative fluxes and energy densities of quantum fields must be subject to 
constraints, at least on average. These constraints, known as quantum inequalities, or QI, provide state-independent 
lower bounds on certain weighted averages of the stress-energy tensor. They may therefore be regarded as the remnants 
of the classical energy conditions after quantisation; thus, for example, the analogue of the WEC is sometimes referred 
to as a quantum weak energy inequality or QWEI. Such bounds typically take the form 

(: T afc ( 7 (T)) U a (r) U 6 (r) :)„ f(r) dr > -0( 7l /) , (1) 

where / > is the weight, or sampling function, 7 is a smooth timelike curve, parametrised by proper time r and 
with four-velocity u a = (d 7 /dr) a , and : T a b : denotes the stress-energy tensor, normal ordered with respect to some 
reference state loq. The significant point is that the bound £2(7, /) is independent of the state u. The class of states 
lo for which the bound holds must also be delineated — all bounds in the literature require (at least) that u> and ujq be 
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Hadamard states, although this can be weakened slightly. Of course the normal ordered energy density differs from 
the renormalised energy density by the (renormalised) energy density of the reference state. Accordingly, Eq. may 
easily be converted into a bound on the renormalised energy density. 

As we will shortly describe, various authors have established QWEI's for the scalar and Dirac fields in different 
circumstances, leading up to general results valid in curved spacetimes. The main aim of the present paper is to 
establish similar QWEI's for the Maxwell and Proca fields in general globally hyperbolic spacetimes. In order to 
do this, we have also made a detailed study of the field theories concerned, particularly in relation to the class of 
Hadamard states. Several new results obtained here may therefore be of more general interest. 

As already mentioned, the earliest work on quantum inequalities is due to Ford, who established a bound on negative 
energy fluxes for scalar fields in 1991 0|. The first QWEI was obtained for the massless, minimally- coupled scalar 
field in Minkowski spacetime by Ford and Roman 0], who established the inequality for the case of an inertial 
worldline 7 and with / taking the Lorentzian form / (t) = To/ {^(j 2 +Tq)), in which To sets the characteristic timescale 
for the averaging. In four dimensions, the resulting bound was 



These results were then extended by Ford and Roman to the massive scalar field in Minkowski space [14l 1 1 5| and 
were further extended by Pfenning and Ford to minimally-coupled scalar fields of arbitrary mass in ultrastatic curved 
spacetimes 0,^3^3- I n these generalisations, it was found that the dominant term of the QWEI has a form similar 
to Eq. (J5J), but with subdominant correction terms due to the curvature and the mass of the field. Pfenning and 
Ford also showed that one could express the bound on the right hand side of the QWEI in terms of derivatives of 
the Euclidean Green's function for the spacetime and developed a short sampling time approximation to the QWEI 
which could be used in spacetimes where it would be too hard to calculate the exact QWEI bound. 

The restriction to the Lorentzian weight to/(t(t 2 + t§ )) was removed by various authors. By making use of the 
conformal properties of field theories in two dimensions, Flanagan pH l20j | has derived optimal quantum inequalities 
for the massless scalar field for arbitrary smooth positive sampling functions, and Vollick has done the same for the 
Dirac field |2l|. QWEI's for the minimally coupled scalar field in static curved spacetimes of any dimension for an 
arbitrary sampling function were established by Fewster in work with Eveson |22| and Teo |23| . 

More recently, techniques drawn from microlocal analysis have been used to considerably generalise previous QWEI's 
and to put them on a mathematically rigorous footing. Fewster [24j used these techniques to derive a QWEI for 
minimally coupled scalar fields in general globally hyperbolic spacetimes (the most general class on which the Klein- 
Gordon equation is well-posed). In this case, averaging takes place along an arbitrary timelike worldline, using any 
weight of the form f(t) = g(t) 2 for g smooth, real- valued and compactly supported. Subsequently, Fewster and Verch 
established similar results for the Dirac and Majorana fields in four-dimensional globally hyperbolic spacetimes |25j . 
Averaging over spacetime volumes has been considered by Heifer j2(| , also in great generality. 

Previous work on QWEIs for spin-one fields has focussed on the electromagnetic field, beginning with the work of 
Ford and Roman 0], who derived a QWEI for the case of Lorentzian sampling along inertial trajectories in Minkowski 
space. More recently, Pfenning [27| has derived a QWEI for the electromagnetic field in static curved s pac etimes with 
arbitrary positive weight sampling functions by using the techniques developed for the scalar field in |22t |23j . Using 
similar techniques, Marecki [23 has derived bounds on the fluctuations of the electric field strength, which are of 
interest in quantum optics. 

In this paper we will adapt the methods of 24] to the Maxwell and Proca fields. This depends crucially on the fact, 
first discovered by Radzikowski 29] for scalar fields, that the class of Hadamard states may be characterised in terms of 
a wave-front set condition on the two-point function. Similar reformulations are known for the Dirac field J30l [3lL l32| 
but there is as yet no full treatment for the Maxwell and Proca fields in the literature. (See, however, [33, yj] for 
[non- microlocal] discussions of Hadamard states for electromagnetism using Faddeev-Popov ghosts.) Our treatment 
of this issue has been influenced to some extent by the forthcoming work of Junker and Lledo I35I . although it has 
been conducted largely independently, leading to some technical differences with their approach |66| . 

The paper is structured as follows. The Maxwell and Proca fields are most elegantly described using differential 
forms; accordingly, we begin in Section [H] with a description of our conventions for differential forms and other 
geometric objects which will appear in this paper. In particular, we delineate the class of globally hyperbolic spacetimes 
to be considered; for technical reasons it is convenient to assume that their Cauchy surfaces are compact and have 
trivial first homology group. This is followed by a brief introduction to microlocal analysis leading to the definition 
of the wave-front set for p-form distributions. 

In Sec. IIIII wc describe the quantisation of the Maxwell and Proca fields in globally hyperbolic spacetimes. We 
adopt an algebraic approach, giving a direct construction of algebras of observables equivalent to those obtained by 
Dimock l3al and Furlani |37j ■ We also define the notion of a Hadamard state for these fields, as a state whose two-point 
function |67| is related in a certain way to a one-form Klein-Gordon bisolution of Hadamard form. Such bisolutions 
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and their microlocal properties have been discussed in detail by Sahlmann and Verch 32] ; we may therefore read off 
the microlocal properties of the Maxwell and Proca two-point functions in Hadamard states. This permits us to apply 
the methods of |2J| to obtain quantum inequalities for these fields in Sec. E] Before this, in Sec. IIVI we specialise to 
the class of ultrastatic spacetimes in order to gain further insight into the abstract definitions of Sec. lIIII In particular, 
we prove the existence of Hadamard states in these spacetimes and use deformation arguments [38( to deduce from 
this the existence of Hadamard states in general globally hyperbolic spacetimes obeying our topological restrictions. 
We also compare our approach with other quantisation schemes, including the Gupta Blculer method. As mentioned 
above, we expect that some of the results obtained here to be of wider interest. 

In Sec. IVII we investigate our QWEI in Minkowski space, and in general ultrastatic spacetimes (modulo the usual 
topological conditions). Simple formulae are obtained for the QWEI bound, which are readily compared with those 
previously obtained for the scalar field. In Minkowski space, the QWEI bound is weaker by a factor of exactly two for 
the Maxwell field (as already noted in and by a factor of three for the Proca field. This is not very surprising and 
simply reflects the number of spin degrees of freedom. In curved spacetime, however, the spin-one and scalar QWEIs 
cannot be related in this fashion. To emphasise this, we explicitly determine the QWEI bound in the Einstein static 
universe, providing a concrete example of our ultrastatic QWEI. Two appendices contain the proofs of some technical 
results required in the body of the paper. 

II. PRELIMINARIES 

Units where H = c = 1 are used throughout. The notation C(f(R n ) denotes the space of smooth, compactly 
supported 68], complex- valued functions on M™. 

A. Geometry and forms 

Spin-one fields on globally hyperbolic spacetimes are most elegantly formulated in terms of differential forms. We 
will follow the conventions of |39j| . which we now briefly summarise for the benefit of the reader. 

Suppose TV is a smooth n-dimcnsional manifold which is connected, boundaryless, orientable, Hausdorff, paracom- 
pact, and equipped with a smooth metric of index s [69j. We denote the space of smooth, complex- valued p- forms on 
N by f2 p (7V); the subspace of compactly supported p- forms will be written Cl p (N). Each p-form may be regarded as 
an antisymmetric covariant p-tensor field and we will occasionally use index notation accordingly. Thus the exterior 
product a A j3 e W +q {N) of a e fl p (Af) and j3 G fl q {Af) is given by 

(p + q)l 

and the exterior derivative d : fl p (JV) — ► tt p+1 (Af) is defined by 

(da) ai ... Qp+1 = (p+ l)V [ai a a2 ... ap+l] , (4) 

where the square brackets denote antisymmetrisation and V a is any connection on Af (d is independent of the choice 
of connection). The Hodge ^-operator is defined uniquely as the map * : ffl(Af) — ► £l n ~ p (AT) such that 

aA*(3= ^a ai ... ap P ai - a ? dvo\ g , (5) 

where dvol g is the positive volume n-form associated with the metric g. In particular, (*) 2 = (— l)p( n ~p)+ s on p- 
forms. By combining the Hodge * and exterior derivative, we may define the coderivative S : £1 P (7V) — > il p ~ 1 (Af) 
(with the convention that d annihilates all zero-forms) by S = (— l) n (p _1 )+ s + 1 * d*, which reduces to S = *d* in a 
four-dimensional Lorentzian spacetime. 

The operations introduced above allow us to define a symmetric pairing (•, •) of p-forms under integration: we set 

(U,V) = I UA*V (6) 

for any W,V £ il p (Af) for which the integral exists. Since d(U A *V) = dU A*V-MA *6V, Stokes' theorem gives 



(SU,V) = (U,dV) 



(7) 
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for smooth (p + 1)- and p-forms U and V whose supports have compact intersection. In this sense the operators d 
and S are dual. 

The Laplace-Beltrami operator is defined as — (dd + Sd), i.e., it is equal to minus the Laplace-de Rahm operator 
d<5 + <5d. Where the manifold JV is a Lorentzian spacctime, with signature H — ■ • • — , the Laplace-Beltrami operator 
is also known as the D'Alembertian and will be denoted by □. We wish to point out that Dimock |36| uses □ to 
denote the Laplace-de Rahm operator; our usage is determined by the convention that □ should have principal part 
g^d^dv, in accordance with typical usage in general relativity. 

Finally, the spaces fig (TV) may be given locally convex topologies [7(j and the corresponding topological duals 
fig (A/ - )' will be called the spaces of p-form distributions on A". (The pairing (•, •) provides a natural embedding of 
fP(A0 in fig(A/)'-) In the case p = 0, we will also use V{M) and V'{M) for fig (A/ - ) and fig (TV)' respectively. The 
exterior derivative and coderivative are defined on these spaces by 

(dW)Cf) - U{SJ) (U G fig(AA)', / e fig +1 (A0) (8) 

and 

(SU)(f) = U(df) (U e fig m, f e fig-'W) (9) 

which extend the definitions given for smooth forms by virtue of the embeddings defined above and the calculation JJJ . 
In a similar way, the Hodge ^-operator may also be extended to a map from p-form distributions to (n — p)-form 
distributions by 

(*«)(/) = (-1)*<"-»>W(*/) (w e ogCAO', / e fig^W) ; (10) 

and it is easily checked that the formula S = (— i)™(p- 1 )+ s + 1 * d* remains true for distributions. 



B. Microlocal analysis and the wave- front set 



Our proof of the quantum weak energy inequality turns on the detailed singularity properties of various distributions 
related to the two-point functions of quantum fields in Hadamard states. The information required is encoded in the 
wave-front set of these distributions, which is defined as follows. (See Ref. |4jj f° r a full presentation.) 

We will define the Fourier transform u of u £ Co°(R n ) using the nonstandard convention 

u(k) = J d n y u( y y k -y, (ii) 

which conforms to the conventions used e.g., in [24j]. The Fourier transform can be extended to scalar distributions 
of compact support by writing u(k) — u(fk) where fk — e lk ' v . Given a cone V C (K. ra \{0}), we will say that u(k) is 
of rapid decrease in V if for each N £ N there exists a real constant Cn such that 

\m\ < (j^jpr Vfcet/ < ( 12 ) 

where |fc| denotes the Euclidean norm of k. 

Smooth compactly supported functions have Fourier transforms which decay faster than any inverse power in the 
whole of K™, but the same is not true for arbitrary distributions u of compact support. A well known example is the 
Dirac <5-function, whose Fourier transform does not decay at infinity in any direction. We define the set of singular 
directions to be the set of all k £ W n \{0} having no conical neighbourhood V in which u is of rapid decrease. 

More detailed information about the singularities of u can be gained by localizing the singular directions. In 
particular, the set of singular directions of u at a point x is defined by 

E x (u)=f|E(xu), (13) 

x 

where the intersection is taken over all smooth compactly supported test functions \ <= C^°(M. n ) with x( x ) 0- The 
wave-front set WF (it) of u is then defined by 



WF (u) = {(x, k) £ E™ x (R"\{0» | k £ E x (u)}. 



(14) 
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The wave-front set can be extended in a natural way to distributions on manifolds. Let (AT, g) be a smooth n- 
dimensional manifold of the type discussed above. Each distribution u in V(Af) has a representative u K 6 T>'(M, n ) in 
each coordinate chart (U, n) (with corresponding coordinates denoted 2 M = k(xY) defined so that 

«i.(A/HsI)=«(/°k) (15) 

holds for every smooth function / compactly supported in k(U) C R™, where \g\ is the determinant of g^ v . The 
wave-front set WF (u) is now a subset of the cotangent bundle T*Af with the property that (x, fc) € WF (u) if and 
only if there is a chart (U, k) about x so that (x, k) € WF (u K ), where (x, k) are the coordinates of (xjz) and u K is 
defined as in l|15|). In fact, it may be shown (see Theorem 8.2.4 and the following discussion in Ref. |40f) that the 
restriction of WF (u) to U is given by 

WF(u)nT*£7 = k*WF(ii k ) 

:= {K*{x,k) : (x,k) £ WF(« K )}, (16) 

where the pull-back k* relates (x, k) € T*A/" to its coordinates by K*(x,k) = (x, k). 

There is a natural extension of the wavefront set to p-form distributions. Let vi (i = 0, ... ,n — 1) be a global 
orthonormal n-bein on N (i.e., the vectors obey g a bvfvj = if 3 and thus rf^vfVj = g ab , where rf 3 is diagonal with 
s entries equal to —1 and the rest equal to +1) and let V % be the dual basis of one-forms: V* = rf' J 'g a bVj- For each 
U G f2g(A/")', we may define the component distributions U^.-.i € T>'(Af) by 

^•••i P (/) = Vhh ■ ■■<„ jMifV A • • • A V Jp ) (17) 
for / G D(Af). Then we define the wave-front set of U to be 

WF(W)= (J WF(W tl ..., p ), (18) 

ii,...,i p 

which may be shown to be independent of the particular choice of n-bein i>i. 



III. QUANTISATION OF THE MAXWELL AND PROCA FIELDS 



A. Classical theory 



We will consider the Maxwell and Proca fields propagating on four-dimensional Lorentzian globally hyperbolic 
spacetimes. Each such spacetime is a pair (A4,g) consisting of a four-dimensional, smooth, real manifold A4, with 
the topological properties listed at the start of Sec. Ill Al together with a smooth Lorentzian metric g with signature 
(+,—,—,—). Global hyperbolicity, which will ensure well-posedness of our field equations, requires that (Ai,g) be 
time-orientable and that M. contain a Cauchy surface E, that is, a smooth spacelike hypersurface intersected precisely 
once by every inextendible causal curve in A4. In fact, one may show that M. is diffeomorphic to K x S. We will 
assume for the most part that £ [also referred to as the spatial section] is compact and that S has trivial first singular 
homology group with real coefficients, i?i(E). This is equivalent to the triviality of the compact support de Rahm 
cohomology group H^(A4) 71]; note that this condition excludes, for example, the 3-torus T 3 as a spatial section. 
The compactness assumption is made for convenience only; triviality of -ffi(E) is inessential for the Proca field, but 
appears to be required in order to establish some of our results for the Maxwell field. 

The classical uncharged spin-1 field of mass M > is a real one-form field A £ ^(M) obeying the Proca equation 

(Sd + M 2 ) A = . (19) 

Applying the coderivative 8 we see that SA = and so any solution to Eq. (|T§|) also satisfies the one-form Klein- 
Gordon equation 

(a + M 2 )A = 0. (20) 
Conversely, any solution to if^UI) satisfying the constraint 

SA = (21) 
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solves the Proca equation. The advantage of the system (|20I21|I is that (|20|l has the hyperbolic principal part q a ^d n d. 



in local coordinates and therefore admits unique fundamental solutions E M : f2j(.M) — * ^{M) respectively [371 li 
such that 

(□ + M 2 ) E M — E M (□ + M 2 ) = 1 (22) 

and supp (E M J) C J ± (supp l 7) for J 6 Vl\(M.). Here, J ± {S) 1 the causal future(+)/past(— ) of S, is defined to 
be the set of points in (Ai,g) that can be reached from the set S C M by a future(+)/past(— ) directed causal 
curve. The operators E M extend to J with suppj7 compact to the past/future and for such t 7, ^4 ± = E M J is 
the unique solution of (□ + M 2 ) A^ = J with supp.4 compact to the past/future. In addition, we introduce the 
advanced-minus-retarded bisolution 

E M = E M - E M , (23) 

which satisfies the homogenous Klein-Gordon equation 

(□ + M 2 ) E M = E M (□ + M 2 ) = . (24) 

We also note that, since d and 6 commute with the Klein-Gordon operator □ + M 2 (or, more precisely, intertwine 
its action on zero- and one-forms) these operators also intertwine the action of E M on zero- and one-forms. 

The fundamental solutions for the one-form Klein-Gordon equation allow us to solve the inhomogeneous Proca 
equation, 

(-*d + M 2 ) A = J, (25) 

with advanced (— ) or retarded (+) boundary conditions. Assuming the existence of a solution and applying the 
coderivative, we find 6A = M~ 2 8J. This allows us to rewrite ()25J) as 

(□ + A/ 2 ) A = (-6d + M 2 ) A - dSA = J- M~ 2 d5J , (26) 



to which A ± = E M {J - M~ 2 &8J) are the unique solutions with support in J^ 1 (supp J~) . Using the property 
^ -E ± l 

J M — M l 



SE M — E M S, these may be shown to be the required solutions to l|25f) . We write 



A± = E M (1 - M- 2 dd) (27) 

for the corresponding solution operators and define 

A M = A M - A+ , (28) 

which will later appear in the commutation relations for the quantised Proca field. We will also use the notation Am 
and E M for the bidistributions defined by A M (f,g) = (f,A M g) and E M (f,g) = (f,E M g) [f,g£ Oj(M)). 
Turning to electromagnetism, the theory is described by a one- form potential A £ ^(Ai) obeying 

SdA = , (29) 

which entails that the field strength T = d.4 obeys Maxwell's equations, 

dJ" = and ST = 0. (30) 

Two potentials A and A' are gauge equivalent, denoted by A ~ A', if A = A' + d\ for some \ £ il°(Ai). Since gauge 
equivalent potentials lead to the same field strength, it is really the gauge equivalence classes [A] of solutions which 
are of physical significance rather than the potentials themselves. We can partially fix the gauge freedom by passing 
to the Lorentz gauge 6 A = 0, in which the Maxwell equations are expressed by the M = case of the system (|20I21II . 
Just as with the Proca equation, the fundamental solutions to ([2H| permit us to solve the inhomogeneous Maxwell 
equation 

-SdA = J, (31) 

where, for consistency, the source J € Qq(A4) is required to obey the current conservation equation 6 J = 0. Defining 
A^ = E~^J, we note that SA^ = SE^J = E^8J = and therefore deduce that A^ is the unique Lorentz gauge 
solution to (|31|l with support in (supp J~) . 
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B. Quantisation: Algebras of Observables 

The canonical quantisation of the Maxwell and Proca fields in globally hyperbolic spacetimes was accomplished by 
Dimock j3(| and Furlani respectively (see also Here, we give a direct construction of suitable algebr as of 

(polynomials in) smeared helds for these theories, isomorphic to those emerging from the constructions of [3|||33. We 
begin with the more straightforward Proca case, for which the algebra of observables will be denoted 21m [M , g) ■ This 
algebra is constructed by first using the set of smooth compactly supported complex valued one-form test functions 
Q,q(M.) to label a set of abstract objects {A(f)\f £ £Iq(A4)} which generate a free unital ^-algebra over C. These 
objects are interpreted as smeared one-form fields 

A(f)<=> (A,f). (32) 

The algebra 21m (A4, g) is defined to be the quotient of 21m by the relations: 

PI. Linearity: A (aft + /3f 2 ) = aA{h) + 0A(f 2 ) for all a,/JeC and ft G ^l(M); 

P2. Hcrmiticity: A(f)* = A(J) for all / G ^(M): 

P3. Field Equations: A ([Sd - M 2 ]f) = for all / G £ll(M); 

P4. CCRs: [A(ft),A(f 2 )\ = -iA M (/i,/ 2 )l for all ft G Q^(M). 

Here, Am is the propagator defined in l)28|l. 

For electromagnetism, quantisation is complicated by gauge freedom. Even in Minkowski space this presents serious 
problems: as shown by Strocchi in the Wightman axiomatic approach, the vector potential cannot exist as an 

operator-valued distribution if it is to transform correctly under the Lorentz group or even display commutativity at 
spacelike separations in a weak sense. Dimock ^3£| circumvented such problems by constructing smeared field operators 
[-A](/) which may be smeared only with co-closed (divergence-free) test functions, i.e., / must satisfy Sf = 0. These 
objects may be interpreted as smeared gauge-equivalence classes of quantum one-form fields: formally, 

[A](f)'='(A,f), (33) 

where A is a representative of the equivalence class [.4]; since Sf = 0, we have (dx, f) = (x>^f) = so t ms 
interpretation is indeed gauge independent. Adapting this idea to our present setting, we start with a set of abstract 
objects {[-4](/) |/ g fll(M) with Sf = 0} labelled only by co-closed test one-forms. As before, we use this set to 
generate a free unital *-algebra over C and define the algebra of observables 2l(A4\ g) to be the quotient of this algebra 
by the relations: 

Ml. Linearity: [A] (afi + 0f 2 ) = a[A](fi) + 0[A](f 2 ) for all a,0 G C and co-closed ft € Slh(M); 

M2. Hermiticity: [A](f)* = A(f) for all co-closed / G ftJ(M); 

M3. Field Equations: [A] (Sdf) = for all / G Slk(M); 

M4. CCRs: [[A]{ft), [A](f 2 )} = -iE {f x ,f 2 )t for all co-closed ft G Sl&M). 

(Note that the one- forms / appearing in axiom M3 need not be co-closed.) 

To see that these algebras are equivalent to those constructed in [3(||33, it suffices to observe that, firstly, the field 
operators of [3^, |32j certainly satisfy the relations above (see Proposition 8 of [3(| and Theorem 3 of [33 ; note that 
there are notational differences) and, secondly, that the algebras we have constructed admit no nontrivial quotients, 
as may be seen by applying the theory of Sect. 7.1 in [4q | [73| | . Accordingly, our algebras are isomorphic to those 

of mil- 

C. Quantisation: Hadamard states 

A state on a *-algebra 21 is a linear functional lo : 21 — > C which is normalised so that w(l) = 1 and has the positivity 
property lo(B*B) > for all B G 21. However, not all states on 21m(A4,<?) and QL(M,g) are of physical relevance — 
many are insufficiently regular to permit the definition of the stress-tensor, for example. We will focus attention on 
the class of Hadamard states, for which the stress-tensor may be defined by point-splitting techniques. The Hadamard 
condition was first stated rigorously for the scalar field by Kay and Wald 0]; more recently, Sahlmann and Verch 
■32] have studied the Hadamard form for wave equations with metric principal part. This does not immediately cover 
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either the Maxwell or Proca equations (which are not hyperbolic). Nonetheless, we may exploit the close relationship 
of these equations to the one- form Klein-Gordon equation to define the notion of Hadamard states for these theories. 

To be more specific, a distribution W £ (f2j(.M) x Qq(A4)) is said to be of Hadamard form if its singular part takes 
a prescribed form in a causal normal neighbourhood Af of some Cauchy surface in M . This requirement is implemented 
by requiring that W — Hk should be C k for each k £N, where Hk is a prescribed sequence of distributions on Af X Af. 
Furthermore, if W is a bisolution to □ + Af 2 (or even a bisolution modulo C°°) then the Hadamard form 'propagates' 
in the sense that W will satisfy the above criterion in any causal normal neighbourhood of any Cauchy surface in M. . 
It follows that the difference W — W' between Hadamard form (□ + M 2 )-bisolutions W and W' is everywhere smooth 
on M. X A4. For our purposes, the second crucial property of a Hadamard bisolution W (first noted in the scalar case 
by Radzikowski |29j |) is that its wave- front set is given explicitly by 

WF (W) =TZ:= k; x', -k') £ t*{M x M) : (x, k) ~ (x', k') and k £ F+j . (34) 

Here, T*(A4 x A4) denotes the cotangent bundle of M. x M. with its zero section excised, and (x, k) ~ [x' , k') if and 
only if k! is the parallel transport of k along a null geodesic connecting x and x , to which k is a cotangent vector at 
x (if x = x' this degenerates to the requirement that k — k' is null). We have also used V x to denote the closed cone 
of future pointing covectors at x. 

We are now in a position to define the notion of a Hadamard state for the Maxwell and Proca fields. Our definitions 
are similar to those employed by Junker and Lledo |35| and (particularly in the Maxwell case) were influenced by 
early versions of their work • 

Proca: A state u> on 21m(-A4><?) is Hadamard if there exists a Hadamard form (□ + M 2 )-bisolution Wm such that 

lo (A(h)A(f 2 )) = W AI (A, (l - M- 2 d<5) f 2 ) (35) 

for all fi £ nl(M). 

Maxwell: A state uj on 2l(.M,(?) is Hadamard if there exists a Hadamard form □-bisolution W such that 

u{[A){h)[A]{h)) = W{f u h) (36) 

for all fi £fl 1 (M) with 5 ft = 0. 

Remarks: a) In neither case is the Hadamard form bisolution uniquely determined by the state. 

b) As stated, these conditions appear global in nature and it is not clear, a priori, that any states satisfying these 
definitions exist in general spacetimes. These concerns will be allayed in Sec. IIVEI where we show that Hadamard 
states exist on general globally hyperbolic spacetimes (subject to our usual topological restrictions on S). A key part 
of this argument is the proof (in Appendix that it suffices for equations (|35|l and i|3t)|) to hold for fi supported 
in a causal normal neighbourhood of a Cauchy surface for them to hold for all fi. In combination with an explicit 
construction of a Hadamard state in ultrastatic spacetimes (in Sees. IIV Bl and HV ( "(I this permits us to apply standard 
deformation arguments |38j to deduce the existence of Hadamard states in general. 

c) Since the two-point function a/ 2 ) (/,/') := u>(A(fx)A(f2)) for a Hadamard state u> of the Proca field is given 
by acting on a Hadamard (□ + Af 2 )-bisolution with a partial differential operator, we may use Eq. I|34() and the 
nonexpansion of the wave-front set under such operators to obtain 

WF(w (2) )C"K (37) 

In the Maxwell case, the two-point function u/ 2 ' (/,/') := w([^4](/i)[^4](/2)) is not a distribution (because it is only 
defined on co-closed test one-forms fi). However, the two-point function of the field strength T{f) := [A](Sf) is a 
bi-distribution on two-forms satisfying 

u(F(f 1 )F(h))=W(6f 1 ,Sf 2 ) (38) 
and therefore has wave-front set contained in 1Z. 

d) The Hadamard condition for electromagnetism was discussed in [33ll34]| in the context of a Faddeev-Popov method. 
Here, the action is modified by the addition of a gauge breaking term and the ghost action, describing the dynamics of 
a complex, anticommuting scalar field. In that context, the electromagnetic potential may be smeared with arbitrary 
one-form test fields and so the corresponding two-point function is a distributional bisolution to the massless Klein- 
Gordon equation (at least in the Feynman gauge). In a Hadamard state (on the algebra generated by the vector 
potential and the ghost fields), this two-point function is required to be of Hadamard form [thus generalising l|36() to 
non-coclosed fi] ; in addition, a Ward identity is required to connect this two-point function to that of the ghost field. 
This ensures that the contributions to the stress-energy tensor arising from the gauge-breaking and ghost actions 
precisely cancel in Hadamard states. 
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IV. HADAMARD STATES IN ULTRASTATIC SPACETIMES 



To provide a concrete example of the foregoing definitions, we now specialise to the class of ultrastatic spacetimcs. 
This will lead to a proof of the existence of Hadamard states for the Proca and Maxwell fields in general globally 
hyperbolic spacetimes as well as providing explicit two-point functions for particular Hadamard states which will be 
used in Sect. IVII We recall that a spacetime (M., g) is said to be ultrastatic if M. = K x £ and g = 1 © —h, where 
(£, h) is a smooth Riemannian manifold. As usual, we will also assume that £ is compact and the homology group 
Hi{Y>) is trivial. We will often denote a point in M. by the pair (t,x), with t G R, x G S. 

Our discussion in this section will proceed as follows. We begin with the construction of a one- form Hadamard 
(□ + M 2 )-bisolution Wm on (A4,g). Next, we construct Fock representations of the Proca and Maxwell fields on 
(M., g) with the property that the Fock vacuum in each case is Hadamard according to the definitions of the preceding 
subsection. In each case, the Wm plays the role of the required Hadamard (□ + M 2 )-bisolution. We also discuss 
the relationship of our procedure to the Gupta-Bleuler approach and the method of quantisation in a fixed gauge. 
Finally, we explain how these results imply the existence of Hadamard states in general spacetimes. Some parts of 
the analysis are relegated to Appendix IB! 



A. Construction of a Hadamard (□ + A/ 2 )-bisolution 



As is well-known, the Klein-Gordon equation in an ultrastatic spacetime is readily reduced to the analysis of an 
elliptic eigenvalue problem. Namely, if £ S ^(Ad) is a static one-form (i.e., independent of the ultrastatic time 
parameter t) then A(t,x) = e~ ZUJt £(x) solves the one-form Klein-Gordon equation if and only if 

K£ = uj 2 ^ (39) 

where the operator K acts on f (x) = £o(s)d< + £e(x) (£o € fi°(£), G fi X ( s )) by 

K£, = ((-A| + M 2 )to) dt + (-A E + Af 2 )£ s , (40) 

with A s and As denoting the scalar and one-form Laplace-Beltrami operators on (£, h). 

To analyse K further, it is helpful to have various inner product spaces in mind. For each p, the space f2 p (£) may 
be endowed with a positive definite inner product 



{u, u)ap(e) = / uA*^v, (41) 
Js 

where *s is the Hodge operator on (S, h). By completing f2 p (£) with respect to the corresponding norm we obtain 
the Hilbcrt space A P (E) of square integrable p-forms on S. The direct sum H = A°(E) A 1 (S) therefore corresponds 
to a Hilbert space of static one-forms, on which the K may be defined as a positive self-adjoint operator (by using 
Eq. H4Ufl to define K on f2°(£) © f2 1 (S) and then forming the Friedrichs extension). Because S is compact, K has 
purely discrete spectrum. 

From a geometrical viewpoint, however, the Hilbert space inner product is not completely natural. It is therefore 
convenient to introduce an indefinite inner product on H by 



((£ I V)) = / U^^&dvohix) = (£o,%)a°(e) ~ (&,t/e)aMe) • (42) 

We will denote the resulting indefinite inner product space (also known as a Krein space) by K. In addition, we 
will say that a set {£j : j G J} of vectors (labelled by the elements of some set J) in K is pseudo-orthonormal if 
((£?' | 'fj')} = =1=1 for each j and | = for all j ^ j'; the set is said to be complete if 

77 = £<&Ui»«Oh»& (43) 

holds for all 77 G K, with the sum converging in the topology induced by the norm of H. We will refer to £ G K as 
'timelike', 'spacelike' or 'null' depending on whether ((£ | £)) is positive, negative or zero, respectively. 
The following result is proved in Appendix iBl 
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Theorem IV. 1. Let £j (j £ J) be a complete, pseudo-orthonormal basis for K such that K£j = o^j^j faj > 0) and 
setA j (t,x)=e- iul i t Z j (x). Then 

w M (h,h) = - Y, ^(&\Zj))(A j ,h)(A],h), (/iSfl 1 ^) (44) 

jeJ:LJj>0 J 

defines a one-form Hadamard (□ + M 2 )-bisolution Wm on (A4,g) which is independent of the particular basis chosen. 
For M > 0, we have 

W M (Jx,h) - W M (f2, fi) = -iE M (fi, / 2 ) (45) 

for all fi £ Qq(.M); in the case M — 0, this holds provided the fi are both co-closed. 

Our task is now to analyse the eigenproblem (|39[1 . We may identify various families of eigenfunctions which are 
more or less convenient for different purposes and which may be combined to form complete pseudo-orthonormal 
bases (again, different bases are useful in different contexts). A typical eigenfunction will be denoted £(A, j), where 
A labels the family to which it belongs and j (which labels eigenfunctions within families) takes values in a labelling 
set J(A). Because K is manifestly positive, the eigenvalues may be expressed as the squares of nonnegative quantities 
u>(X, j): thus K£(\, j) — u>(X, j) 2 £(A, j). The corresponding Klein-Gordon solution is 

A(\,j)(t,x) = e-^ x <» t tt\j)(jr). (46) 

We begin by identifying three eigenfunction families which together form a simultaneously K-pseudo-orthonormal and 
H-orthonormal basis. The analysis of Eq. (|39[1 is considerably simplified by the fact that it decouples into the two 
equations 

(-A! + M 2 Ko = u 2 Zo (47) 
(-A s +Af 2 )£ s = u 2 &- (48) 

Accordingly, choosing ipj (j G J(S)) to label a complete orthonormal basis of (— A|. + M 2 )-eigenfunctions for A°(£) 
with corresponding eigenvalues uj(S,j) 2 , we may identify a family of "scalar" K- eigenfunctions given by 

^S,j) = Vj dt (jeJ(S)). (49) 

(By elliptic regularity, each ipj is in fact smooth.) Clearly the £,(S,j) are pseudo-orthonormal and 'timelike' in K, as 
well as being H-orthonormal. For future reference, we note that there is a unique spatially constant eigenfunction 
volft(S) _1 / 2 dt with eigenvalue M 2 . Writing $j(t, x) = e^^^-^ipj^), the corresponding one-form positive frequency 
modes are 

A(SJ) = $ j dt. (50) 

The remaining pseudo-orthogonal modes must have vanishing dt component; that is, they must lie in the 'spacelike' 
subspace A 1 (S) of H. On this subspace, the inner products of H and K differ by an overall sign only, so K-pseudo- 
orthonormality and H-orthonormality are again identical. Since Hi(T,) is trivial, the Hodge decomposition (see, e.g., 
Prop. 11.7 in gives A X (S) = ds^°(S) © <5j]il 2 (I]), where the decomposition is orthogonal with respect to (•, -)ai 
and hence ((■ | ■)), and the bars denote closure. As is easily verified, the Laplacian As is block diagonal with respect to 
this decomposition, thus enabling us to seek eigenfunctions within each subspace in turn. If the restriction on i?i(S) 
were removed, there would be a third subspace, consisting of harmonic forms (i.e., the kernel of As). In the Proca 
case, the resulting modes lead to additional terms in various expansions given below and in the quantum inequality, 
but do not significantly alter the formalism. By contrast, the incorporation of the harmonic modes in the Maxwell 
case raises nontrivial issues, to which we hope to return elsewhere. 

Now, the "longitudinal" subspace ds^°(S) is clearly spanned by the set of non- vanishing vectors of the form 
d£(<5,j) (j £ ^("SO); owing to the relation A^ds = dsAf., these must also be eigenvectors for —As + M 2 with 
eigenvalue uj{L,j) 2 = u>(S,j) 2 . As the only vanishing vector of this form is obtained from the spatially constant 
mode, the appropriate labelling set is J(L) — {j E J{S) : uj(S,j) > M}. Furthermore, the calculation 

((dz<Pj I dx<Pj>)) = -(ds</Ji,dsyy)Ai(£) = -(<pj, 5 s d s ^/) A1(s) = -(u(S, j) 2 - M 2 )S jr (51) 

shows that the appropriately normalised longitudinal eigenfunctions are 



i{L,j) = (u(S,j) 2 - Af 2 )- 1 / 2 d s ^- , (j e J(L)) . 



(52) 
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The corresponding Klein-Gordon modes may also be expressed in the form 

+ iu)(L,j)$j dt 



ML, 3) 



^uj{L i3 ) 2 -M 2 



(53) 



All remaining pseudo-orthonormal modes must lie in the coexact subspace S^fl 2 (Y.). We will refer to these as the 
transverse modes £(T,j) with labelling set J(T) and eigenfrequencies u)(T,j); they are necessarily 'spacelike'. In 
general, the ui(T, j)'s will be distinct from the u>(S,j)'s (even in situations of high symmetry such as the Einstein 
spacetime — see Sect. IVlUj) . 



Applying Theorem lIV.ll we obtain a Hadamard bisolution 

S (A) 



w M {h,h) = - E E 



\e{S,L,T}j£J(\) 



(A(X,j),h)(A(X,j),h), 



(54) 



where s(S) = 1, s(L) = s(T) = -1. 



The Proca field 



We now construct a state on QlM(M.,g) whose two-point function is related to Wm by Eq. (|35|l and is therefore 
Hadamard. We begin by finding an expression for Wm(Ji, (1 — M~ 2 d5)f2). To this end, we first note that the 
transverse modes £{T,j), being coexact, are necessarily coclosed: S^(T,j) — 0. It follows that SA(T,j) — 0, so 
transverse modes therefore solve the Proca equation. The same cannot be said of the scalar and longitudinal modes, 
for which 



6A(S,j) = iu(S,j)$ j 



for all j £ J(S) , and 



8A(L,j)=-^J{Sjf 
for all j £ J(L). It is convenient to change basis, introducing 

1 



£(P,j) - - [u(S,j)t(L,3) - Wu(S,j) 2 - M*&S,j) 



= [V^(S,J) 2 ~ M*£(L,j) - iLj(S,j)Z(S,j) 



(55) 
(56) 

(57) 
(58) 



with labelling sets J(P) = J(G) = J(L). The two sets are pseudo-orthonormal, with £(P,j) and £(G,j) being 
'spacelike' and 'timelike' respectively. By construction, the "scalar Proca" modes A(P,j)(t,x) = e~ %ul ^ p ^ t ^{P, j){x) 
are coclosed and therefore obey the Proca equation. We also note the expressions 



A(PJ) 



1 

M 



u>(S,j)A(L,j) - iy/u(S,j) 2 -M 2 A(S,j) 
iM 



On the other hand, the corresponding "gradient" modes A(G,j) have 

A(G,j) = M~ 1 d* j 

and therefore have nonvanishing coderivative 

SA(G,j) = M$ j . 



M 2 



(59) 



(60) 



(61) 



From Theorem lIV.ll we know that Wm is basis independent and may therefore be expressed in terms of the A{P, j), 
A(G,j), A(T,j) and the unique spatially constant mode. Of these, the scalar Proca and transverse modes are left 
invariant by the operator (1 — M~ 2 d5), while the spatially constant mode and the gradient modes are annihilated. 
Accordingly, with Wm expressed in the new basis, we find 



W M {{1 - M- 2 dS)hJ 2 ) = W M (fi, (1 - M- 2 d<5)/ 2 ) = J2 E 



1 



Xe{T,p}jeJW 



2w(A,j 



■(A(X,j),fi) (A(X,j),f2) 



(62) 
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As a consequence of the first equality and Eq. (|45|l we see that 

W M (fi, (1 - M- 2 dS)f 2 ) - W M (h, (1 - M- 2 d8)h) = -iE M (.fu (1 - Af- 2 d^)/ 2 ) = -iA M (/i, h) (63) 
which in turn leads to the expansion 

A M (/i,/ 2 )= Y, E ^^A(A{\])Ji)(J^)j2)-{A{\ ] )J 2 ){M^ 1 h)) ■ (64) 

To complete the analysis, we define smeared fields 

A W = E E Ar 4^-r (W,j),/)a(A, j) + (A(XJ)J)a(X,jr) (65) 

on the Fock space Fm generated by operators a(X,j) obeying [a(X,j),a(X',j')*] — 5\\>8jji1L. It is now easily verified 
that the A(f) provide a representation of axioms P1-P4 on (a dense domain in) Fm- Thus the Fock vacuum |0) 
defines a state on 21m (-M, g) with two-point function 

(0\A(h)A(fr)\0) = W M (fu (1 - M- 2 dS)f 2 ) . (66) 

Accordingly, this state is Hadamard. Acting on the vacuum with smeared fields, we obtain a dense set of Hadamard 
states in Fm- 



C. The Maxwell field 



We now repeat the above analysis for the Maxwell field. In this instance, we begin by finding an expression for 
Wo(/i, f 2 ) for co-closed /,-. Since the A(P,j) and A(G, j) diverge as M — > 0, we revert to the scalar, longitudinal and 
transverse modes of Sect. llV Al Now, from Eqs. (p)0|) and it is simple to show (A(L,j),F) — i(A(S,j),F) for all 
co-closed F. This has the effect that the contributions to Wo(fi, f 2 ) from the scalar and longitudinal modes cancel 
and so only has contributions from the transverse modes: 

W (h,f 2 )= Y ^J^{A{T,j)J{){MTj),f 2 ), (67) 

for any co-closed fi. Defining 

[A](f)= Y Ar 4^ ((A(T,j),f) a(T,j) + (AtTjj, f) a(T,j)*) (68) 

on the Fock space Fo generated by operators a(T,j) obeying [a(T,j), a(T, j')*] — Sjj'l, we may again verify that the 
[^4] (/) provide a representation of axioms M1-M4 on (a dense domain in) F . Thus the Fock vacuum |0) defines a 
state on 2l(.M,<?) with two-point function 

(0|[^](/i)[^](/ 2 )|0) = WbCfi,/2) (69) 

for all co-closed fi and is therefore Hadamard, as are states obtained from it by acting with polynomials in smeared 
fields. 



D. Comparison with other quantisation methods 



It is instructive to compare our approach to two more familiar methods of quantisation, in particular, the Gupta 
Bleuler formalism [fifl l51| and (for electromagnetism) quantisation in a fix ed g auge. A detailed presentation of these 
methods for electromagnetic fields in static spacetimes can be found in [27fl. We will confine our remarks to the 
Maxwell case, but analogous comments apply to the Proca field. 

As is well-known, the Gupta-Bleuler procedure starts by adding a gauge breaking term to the Maxwell Lagrangian. 
By choosing the coefficient of this term appropriately, the field equations become precisely those of the massless 
Klein-Gordon equation. (This is sometimes known as the 'Feynman gauge'.) This equation is then quantised using 



13 



an indefinite inner product space (effectively the symmetric Fock space F(K) over our Krein space K). The physical 
Hilbcrt space is then selected by the requirement that they be annihilated by the positive frequency part of the 
divergence of the resulting field operators. 

In our approach the various ingredients of the Gupta-Bleuler formalism appear in a different way. We start from 
the Hadamard bisolution Wo, which is in fact the two-point function of the static vacuum of the Gupta-Bleuler 
theory in the indefinite space F(K). In our approach, however, the restriction to the physical state space occurs at 
the 'one-particle level': the Hadamard condition and the restriction to co-closed test one-forms pick out a preferred 
Hilbert subspace of K, namely, the subspace H p h ys spanned by the transverse modes £(T,j). The quantum fields are 
then given immediately as operators on the symmetric Fock space Fo = F(H p h ys ) over this Hilbert space. The key 
advantage of our approach over the Gupta-Bleuler method is, of course, that it is not tied to a particular notion of 
'positive frequency' and is applicable in general globally hyperbolic spacetimes. 

An alternative to the Gupta-Bleuler approach is to fix a gauge from the outset. Employing the Coulomb gauge |27| . 
we would seek modes Aj obeying SAj = with vanishing dt component. These modes are easily seen to constitute 
the transverse family A(T, j) and we essentially recover the field operators of Eq. (|68|l . but without the restriction to 
co-closed fi. By contrast, no gauge is ever chosen in our approach and the transverse modes emerge naturally from 
the Hadamard condition in combination with the restriction to co-closed fi. 



E. Existence of Hadamard states in globally hyperbolic spacetimes 



Once the existence of Hadamard states has been established in ultrastatic spacetime, arguments due to Fulling, 
Narcowich and Wald ma y be used to deduce their existence in general globally hyperbolic spacetimes obeying 
our usual topological restrictions. We illustrate this for the Proca field, indicating the slight differences required in 
the Maxwell case. 

Let {M,g) and (A4\g') be globally hyperbolic spacetimes and suppose that there is an isometry ip : N — > N' 
between causal normal neighbourhoods Af and Af' of Cauchy surfaces in M and M' respectively. Then any state w on 
%m(A4, g) determines a state ui' on $Im{A4', g') as follows: given any /i, . . . , /„ € fll(M') choose fx, . . . , /„ € fig (A/ 7 ) 
such that fk — fk S (— Sd + M 2 )fll(A4'). (The existence of such fk follows from Proposition IA.3t aV) Now define 

u'{A\h) ■ ■■A'Un)) = u{A{rJx) ■ ■■A{rJn)) , (70) 

where ip* denotes the pull-back, ip* f = f o ip and we have used A'(f) to denote field operators in ^.m(A4', g'). 
It is easily verified that u>' is a state on 21m (A4' , g 1 ). Moreover, if u> is Hadamard then the isometry ensures that 
Lj'(A(f 1 )A(f 2 )) = W(fi, (1- M~ 1 dS)f 2 ) for all f t G ^(A/ 7 ), where W is a (□ + M 2 )-bisolution Hadamard bisolution 
on Af x Af. It then follows from Theorem I A . 1 1 that uj' is Hadamard. 

Given a globally hyperbolic spacetime (A4 : g), we may now employ a construction described in |38| to obtain a 
globally hyperbolic spacetimes (M',g') and (M",g") so that (i) (M",g") is ultrastatic; (ii) M' contains Cauchy 
surfaces E^ and Ef, with causal normal neighbourhoods M[ and M 2 so that N[ (respectively, Af 2 ) is isometric to a 
causal normal neighbourhood of a Cauchy surface in (A4,g) (respectively, (A4",g")). As all the Cauchy surfaces 
involved are therefore homeomorphic, (A4',g') and (A4",g") will obey our topological restrictions if (A4,g) does. 
Starting with a Hadamard state on the ultrastatic spacetime (A4",g") we may induce Hadamard states on (A4',g') 
and hence on (M , g) . Thus 21 m (M , g) admits Hadamard states. 

An analogous argument applies in the Maxwell case. The only differences are that all test functions must now 
be co-closed, we use part (b) of Proposition IA.3I instead of part (a), and use the appropriate form of the Hadamard 
condition. 



V. A QUANTUM WEAK ENERGY INEQUALITY 

We now proceed to prove our quantum inequality. Let r i— > 7(7") be a smooth timelike curve in M. parametrised by its 
proper time t£K, and let T be a tubular neighbourhood of 7. Inside T let {v? \i = 0, 1, 2, 3} be an orthonormal frame 
obeying g pa — rf^v^vj and with the prop erty that the restriction of Vq to 7 is the four-velocity u a (T) — (^(t) / dr) 17 
of the curve: Vq (7(1")) = u CT (r). (See [25| for an explicit construction of such a frame.) 

The first step is to construct the quantised energy density measured along the curve. The starting point is the 
classical stress-energy tensor, which for the Proca and Maxwell fields, takes the form 

-^L*i/ — ~^g^u^~ p<jJ~ p J~ (tpF ~\~ M I A^Av —g ' ^ii/ApA^ J , (71) 
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with M = in the Maxwell case. The energy density along the curve is given by T M ^(7(r))w^ (t)u u (t); we may extend 
this quantity off the curve by defining T{x) = T^ u (x)vq (x)vq (x). A little manipulation shows that 



3 3 3 



r (-) = z E E K«^) 2 + \ m2 E «a>) 2 • ( 72 ) 



i=0 j=0 i=0 



We may rewrite T in terms of forms with the aid of the one-form basis {V 1 : i = 0, 1, 2, 3} dual to the Uj, given by 
V£ = Writing = A V, we find 



3 3 .3 



T w 4EEK fA * wij )] 2 + ^ E [* M A *^)] 2 . (73) 

i=0 j=0 i=0 

which is clearly the restriction to the diagonal x' = x of 

T (*< = iEE [* L [* ^ A *^)] L< + \ M2 E [* M A *^)] I* [* M A *^)] I • ( 74 ) 



i=0 j=0 



The function T(x, x') provides a point-split classical energy density, related to the true energy density on the curve 7 
by T( 7 (r)) = T( 7 (t), 7 (t)). Moreover, for /, /' E V(M) we see that 



3 3 



/ T(x, x')f{x)f(x')dvo\ g {x) dvol g (x') = - W^f){T, W ij f) + \m 2 ^{A, V l f){A, V*f) 

J MxM 4 4 =0j=0 1 i=0 

3 3 3 

7 E EM. &{W ij f)) M, 5{W ij f')) + ~M 2 J2(A, V l j){A, V'f) 



i=0 j=0 i=0 



(75) 



The advantage of this reformulation is that it is easily quantised: given a Hadamard state co on 21m (-M, 9) for the 
Proca field, or Ql(A4,g) for the Maxwell field we simplify replace occurrences of (A,F)(A,F') with the two-point 
function lo^(F,F') thus obtaining a scalar bi-distribution (T) w e 2?'(.M x A4) given by 



33 3 
(T)MJ') = JEE^ (2) (*(W* j 7), S{W ij f)) +Im 2 <TujW (V'f,V l f) , (76) 

i=0 j=0 i=0 

(Note that in the Maxwell case, the above expression reduces to its first term, in which the arguments of the two-point 
function are co-closed.) 

Now suppose that a reference Hadamard state u>o is specified. Since the difference between Hadamard form (n+M 2 )- 
bisolutions is smooth, 

(: T :) u = (T) u (T) Uo (77) 

is easily seen to be a smooth function on M x A4. Thus we may 'unsplit points' to define the normal ordered energy 
density along 7 by 

(:p:) u (r) = <:T:) w ( 7 (r),7(r)). (78) 

We will need three further observations. First, (: T :} ul is symmetric as a consequence of the commutator axiom 
P4/M4. Second, {T) u is a distribution of positive type, i.e., (T) u (f,f) > for all / e V(A4). This is ultimately a 
consequence of the positivity property uj(B*B) > and the hermiticity axiom P2/M2. Third, the wave-front set of 
(T) u obeys 

WF((T) U )CTZ. (79) 

In the Proca case, this is a straightforward consequence of the non-expansion of the wave-front set under partial 
differential operators. In the electromagnetic case, the key point is that (F,G) 1— > lj( 2 \5F, SG) is a distribution (the 
two-point function of the field strength) whose wave-front set is known to be contained in 1Z. 
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Now consider the the pull-back <^*(T) W of the bi-distribution (T) w induced by the smooth map (p(r,r') = 
(7(7"), 7(7"')). This quantity is formally the unrenormalised energy density, with points split along 7. We now 
claim that ip* (T) w is a well defined distribution of positive type. To verify this, we first note that t <p'(r,T') : 
T* 7 ( r ) J f T n\(M x M) — > K is the linear map 

V(t,t') : ^ {uP{r)k p y {r')k' a ,) (80) 

from which it follows that p has the following set of normals 

N<p = Ul{T),k;i(T'),k')eT*(M xM) :u p {T)k p = u a ' {T r )k' a , = o} . (81) 

Second, we note that all the covectors appearing in WF {{T) u ) are null, and therefore cannot annihilate any nonzero 
timclike vector. Accordingly, the intersection WF ((T)^) n N v is empty and the pull-back tp*(T)u exists by Theo- 
rem 2.5.11' in 52] (in which the set of normals is also defined). Theorem 2.2 of [24j guarantees that it inherits the 
positive type property. 

Theorem 2.5.11' in [12 also asserts that the wave-front set of obeys WF {ip*(T) u ) C p*WF ((1%). Thus, 

we have (r, C; t', -(') e f'WF ((T)J only if 

(C-C) - (V'(r,r')) (*,-*') = K(r)fc p ,-^'(r')fc',). (82) 

Since the vectors u ct (t), u ct (r') and the covectors k a , k' a , are all future pointing their contractions will always be 
positive. It will therefore be the case that (, £' > 0. 

Summarising, (y9*(T) w is a well defined distribution in V(R 2 ), which is of positive type and has a wave- front set 
obeying 

WF (<p*{T) u ) C {(r, (; t', -C')l C, C > 0} . (83) 

We are now in a position to state and prove our quantum weak energy inequality. The proof is in fact essentially 
identical to that given for the scalar field in |2J|; it is given here for completeness. 



Theorem V.l. Letuj andujo be Hadamard states on 21m (A4, g) for the Proca field, or 2l(.M , g) for the electromagnetic 
field. Define the point-split normal ordered energy density (: T :) w relative to luo by Eq. |77| ) and the normal ordered 
energy density (: p :) w (r) by Eq. \7$jj . Then the quantum inequality 

dr(g(r)f(:p:Ur)>- — [(g ® g)<p*(T)„ ] A (-a, a) (84) 
Jq n 

holds and the right-hand side is finite for all real-valued g € Cq°(M). 

Proof. Letting g G Cq°(]R) be a real valued function of proper time r, we have 



dr(g(r)r(:p:)Ur) = / dr (<?(t))V<: T ^(r, r) 

dr f dr' g(T)g(T>)S(T~T')ip*(:T:) u (T,T>) 

dr J dr' g^ a (T)g a (r')p* (: T :) u) (t,t / ) 

/>OC j 

= / — ¥>*<:T:) w (^® 5a ) (85) 
Jo n 

where g a {r) = g(r)e laT and we used the symmetry property T (8 .g Q ) = <ys*(: T :) w (g Q ® in the last 

step. Using the definition of (: T :) w and the fact that ip*(T) u is of positive type we have 

dr ( 9 (t)) 2 p*(:T:Ut,t)>- — ^(T)^ <g> fla ). (86) 



da 
2tt 
do 
2^ 
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Setting ei a , a ') ( T i T ') — ej( aT+a T \ we note that the integrand in the above expression may be rewritten as 

f*(T)^ (g^(g) g a ) = <p*{T) U0 ((g® g)e(_ aiO1 y) 
= [{9® g)^P*{T)uj ] (e(- a , a )) 

= [(g®g)v*(T) U0 ] A (-a,a) (87) 

by our definition for the Fourier transform. Thus we obtain Eq. I|84|l . 

The fact that the right hand side of Eq. I|84l) is convergent follows from an analysis of the wave- front set. From the 
definition of the wave- front set, it is obvious that 

WF ((A ® h)<p*(T) Uo ) C WF ((p*(T) Uo ) (88) 

for any /, e Cq°(K). Thus the singular directions for (/i (g) f2) ( P*(T} LJo are contained in {(C,-C')IC)C > 0}- In 
consequence, the Fourier transform [(g <g> g)ip* (T) Wo ] A (—a, a) decays rapidly as a — > +oo. The integral on the right- 
hand side of Eq. (|84|l therefore converges thereby providing a finite bound for all g £ (R) . □ 



VI. EXAMPLES 



A. Minkowski Spacetime 



Although Minkowksi space lies outside the discussion of Sec. IIIII because its Cauchy surfaces are noncompact, 
the arguments presented in Sec. [V] apply equally well to the two-point functions arising from standard Minkowski 
quantisations of the Maxwell and Proca fields based on the Poincare invariant vacua. As we will see, these vacua 
are Hadamard; furthermore the topology of the Cauchy surface played no role in the proof of our QWEIs. We refer 
the reader to .53J (and references therein) for a careful discussion of the quantisation of constrained systems in the 
context of algebraic field theory in Minkowski space. 

We illustrate our QWEIs for the case of an inertial worldline 7, using the usual Poincare invariant vacua as 
reference states. Boosting to the rest frame of 7, we adopt coordinates in which j(t) = (r, xq, yo, zq) for some fixed 
(xo, Ho, z o) £ R 3 - We also adopt the orthonormal frame = (d/dx 1 )^ for i = 0, . . . , 3. 

Let us begin by recalling that the vacuum two-point function for the scalar Klein-Gordon equation is given by 

where k j = [(J,k) and uj = (|fc| 2 + Af 2 ) 1 / 2 . This is of course a Hadamard (□ + A/ 2 )-bisolution, and it may be used 
to define a one- form Hadamard (□ + M 2 )-bisolution by 

Wu(f,f') = -V ij W^(f i ,f j ). (90) 

(The overall sign is determined by the requirement that — anticipating the commutation relations P4/M4 — the anti- 

symmetric part W M (f, f) - W M (f, /) should be equal to -iE M (f, /') = -^'^ ) (/ l , /,)■) 
Turning to the Proca field, the Poincare invariant vacuum two-point function |54| 

-fC/,/') - - [ ( V ij - %{-k)l(k) (91) 



/ R 3 (2tt) 3 2lu V M 2 
is easily seen to be of the form 

4 2) (/,/') = W M (f, (1 - M- 2 A8)f) (92) 
and is therefore Hadamard. We may also write the vacuum two-point function in terms of its kernel 

4 a *W) = -/ Si^-^V- 1 "'''- (93) 



/ R 3 (2tt) 3 2lu \ ' M 2 

It is not a difficult calculation to show that the point-split vacuum energy density along the worldline 7 is 

^(TU(t,t>) = \ / R3 70^" MT - r,) = 7^)3 jf^^Me^' , (94) 
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where, in the last step, we have changed to spherical polar coordinates in the momentum integration, performed the 
angular integrals, and written oj(k) = \/ k 2 + M 2 . This expression can be used to evaluate the right-hand side of 
Eq. giving 

dr (: p :Ur) 9 (r) 2 > - / — [(g <g> g)<p* (X% ] A (-a, a) 
Jo * 

da / dn k 2 lj(k) \g (a + uj(k))\ ■ 



(27T) 3 Jo 

Before commenting on the above expression, we wish to perform a similar analysis for the electromagnetic field. 
The vacuum two-point function for electromagnetism in the Coulomb gauge is given by 

4 2)ij (^') = -/ 7^3^^e-*(— ')\ (95) 

where to = and 

h H = { il lJ +k l W/u 2 fori^O^j 
1 otherwise. 

It is straightforward to check that, for co-closed test one-forms /, /' [i.e., k l fi(k) — 0] we have ui (/, /') = Wo(f, /') 
and that loq is therefore Hadamard. An identical computation to the above now reveals that 

^(TU(t,t<) = / ^ue-M'-O = /°° d KK 2 w(K)e --(«)(r-r') ; (97) 

and hence 



(2tt)3 (2tt)3 7 



dr( :j0 :) w (r) 5 (r) 2 >-— - / da/ d« \g (a + . (98) 

\ Z1T ) Jo Jo 

From this point on, the Proca and electromagnetic fields can be treated simultaneously as the two bounds take the 
common form (setting M — for electromagnetism) 

dr{: P -) u {T)g{T) 2 >-sI M (g), (99) 

where s is the number of spin degrees of freedom for the field, i.e. s — 2 for Maxwell and s = 3 for Proca, and 

2 

(2*) 



M<7) = 7^3 / da I dKK 2 u{n)\g{a + Lo{K))\ 2 . (100) 



This integral can be further simplified by making an additional change of variables, 



u = a + w(k), 

v = uj{k), (101) 



whereupon the above integral can be rewritten as 

2 



M JM 



/A/(ff)= (2^' du\g(u)\ 2 I <WvV-M 2 . (102) 



In fact, the scalar field also satisfies Eq. (|99|) with s = 1 |22|. Thus, in Minkowski space, the QWEI bounds for 
scalar, Proca and Maxwell fields differ only by a factor equal to the number of spin states (as already noted for the 
electromagnetic field in 27]). However, in more general curved spacetimes the scalar and one form fields can have 
different eigenfrequency spectra leading to very different quantum weak energy inequalities for each field. 
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B. Ultrastatic Spacetimes 

We now turn to the class of examples in which the spacetime is ultrastatic (with compact Cauchy surface obeying 
our usual requirements) and the energy density is sampled along a static worldline, using the ultrastatic vacua of 
Sect. II VI as reference states. In this case, the quantum inequality Eq. (|84|) takes a remarkably simple form. Recall 
that, for both the Proca and the Maxwell fields, the two-point function of the vacuum state can be written compactly 
as 

J 2 \F,F')= J2 E ^- ] (A(\j),F)(A(XJ),F') (103) 

A physical j6J(A) ^ 

where A ranges only over the physical families of modes, i.e. A £ {T, P} for Proca, but A = T only for Maxwell. In 
the Maxwell case, we restrict attention to coclosed test one- forms F and F' . Using the mode expansion l|103|) for the 
two point function, and recalling that T = d„4, it is not difficult to show that [74( 

l(g®g)v*(T)X(~a : a)= J ' dr j ' dr' g(r) g(r') e~ ia ^ ]T £ N(X, j)( 7 (r), 7 (r')) , (104) 

A physical jgJ(A) 

where 

i fi 3 3 r n 

N(\,j)(x,x') = —— -EE[V(AJ>MW ^ pCT (A,j>X W) 

+ 1 -m 2 [M^JK] (*) [M^TH] (*')} ■ (105) 
p=o J 

Now if we choose a static worldline 7 (t) = (r, Xo), then 

^(A,j)( 7 (r), 7 (r')) =e"- ^>^'> £(A, j) , (106) 

where 

{ 3 3 3 ^ 

i EE I (7(ro))| 2 + \M 2 I [MMW] (7(ro))| 2 (107) 

is the classical energy density per mode at the spatial position of the observer's worldline. (Note that the leading 
factor of (2w(A, j)) -1 appears because our modes were normalised using the L 2 -inner product on E, rather than the 
symplectic inner product.) The most notable point to be made here is that T(A, j) is r independent because the e~ lUxT 
dependence in the modes is removed by taking the magnitude of the complex functions. Thus eq. (|1U4|) becomes 

[(g^g) lf *(T)^(-a,a) = £ £ \g(a + w(A, j))\ 2 T(A, j) , (108) 

A physicalj'G J(A) 

and the quantum weak energy inequality for static observers in ultra-static spacetimes becomes 

dr(g(T)) 2 (:p:Ur)>- — E E l?(« + ^(A, j))| 2 T(A, j). (109) 

J J ° ^ A physical j£J(A) 

As we will see, the above form turns out to be calculationally convenient in some settings. However, we also note that 
the quantum weak energy inequality may be written 



dr(g(r)) 2 (: p :)„(r) > - / du\g(u)\ 2 Q{u) (110) 



where 



g(«) = - E E s ( A 'j) ( m ) 

A physical 

s.t. k>(A,j)<ii 
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measures the maximum e nerg y density available if no mode with frequency greater than u is excited. The analysis 
carried out in Sect. 5 of |24j for the scalar field may be adapted to show that Q(u) is positive and polynomially 
bounded as u — > oo. Several other systems obey quantum inequalities of this form, which motivated the study of 
related conditions for general quantum dynamical systems in |55|. In the terminology of that paper, the ultrastatic 
vacuum states constructed here fulfill a limiting static Q WEI. We expect that the analysis of [53 would apply equally 
well to the Proca and Maxwell fields, thereby establishing links between quantum inequalities and the thermodynamic 
properties of these fields. 



C. Static Einstein Universe 



A nice example of the ultrastatic quantum inequality is provided by the static Einstein universe Rx S 3 with line 
element 

ds 2 = dt 2 - a 2 [d X 2 + sin 2 X (dO 2 + sin 2 9dip 2 )] . (112) 

Here a is the radius of the universe and (x, 0, ip) are spherical polar coordinates on the unit three-sphere. Various 
authors have studied the electromagnetic mode functions in the Einstein universe |56l 1571 158| and it turns out that 
essentially the same ansatz can be used to obtain the modes of the Proca field. 

We know from Section ITvl that the scalar Proca modes are completely specified by finding mode solutions to the 
massive, minimally-coupled scalar wave equation (o s + A/ 2 ) $ = 0. A complete set of suitably normalised positive 
frequency mode solutions in the Einstein universe is given by |59l l60| 



t>ni m (x) = a~ 3 / 2 n n i( x ) Y lm (6, <p) e- to »*, (113) 

where Yi m (9, tp) are the standard spherical harmonics on a two-sphere |6l| . H n i(x) can be written in terms of Gegen- 
bauer polynomials C^{x) j^] as 

H ni (x) = // . === - sin< X C^(cos X ), (114) 

yjiryn + l + ly. 

and the frequency of the modes is 



n(n + 2) , 

' + M 2 . 115 

Here, the primary quantum number n ranges over the non-negative integers n = 0,1,2,.... For a given n there 
are (n + l) 2 harmonic states with the same frequency labelled by the quantum numbers, I = 0, 1, . . . , n and m = 
— /, — I + 1, . . . , 0, . . . , I — 1, 1. However, the case n = I = m = corresponds to the spatially constant mode; thus the 
labelling set for the scalar Proca modes given by Eq. (|59l) is n = 1, 2, . . ., with / and m as before. 

On the other hand, the family of transverse modes breaks up into two subfamilies |58(: the magnetic J-pole modes 
taking the form 

A(M, n, I, to) = *d (^ n i m dt A d*) , (116) 

and the electric J-pole modes, 

A(E, n, I, to) = *d [-A(M, n, I, to) A dtj , (117) 
where the scalar functions ^ r dm obey 

4^d x + Af 2 ) V nlm = 0. (118) 
a 2 - sm x J 

The positive frequency mode solutions to this partial differential equation are 

* nlm (x) = V2o7lV rd (x)Y lm {e,ip)e- l ^ t , (119) 

where 



VAX) = ,/ ^ / ) ; [ sin'+7e i i - 1 (cosx) (120) 

y / Trl(l + l)(n + l)\ 
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and the frequency of the modes is given by 



UJ n = \ — + M 2 . 



(121) 



Here the quantum numbers range over a different set of allowed values. The primary quantum number n ranges over 

the integers n = 2,3 For a given n there are n 2 — 1 harmonic: states with the same frequency labelled by the 

quantum numbers, I = 1, 2, . . . , n— 1 and to = — I, — 1 + 1, . . . , 0, . . . , I— 1, 1. We stress that both transverse sub- families 
have an energy spectrum given by Eq. (|121f> . which differs from that of the scalar Proca modes. There is actually a 
more convenient basis of transverse mode solutions given by 



A(±, n, I, m) = —= \A(M, n, I, to) ± A(E, n, Z, to) 
v2 



(122) 



which are simultaneous eigenfunctions of (5s *s and — As + M 2 (which commute). The A = M or E families of modes 
are, roughly speaking, the linearly polarized solutions of the field while the A = + or — modes are the circularly 
polarized solutions. 

For the interested reader, we also state the component form of the three families of physical modes. The two 
transverse polarizations of the field are 



A^{M,n,l,m) 



1 



0, 0, 



sm( 



-<%,, sindde I >F 



(123) 



and 



Aa(E, n, I, m) 



0. 



1(1+1) 



sin x 



- nlm ■ 



(124) 



while the scalar Proca modes are 

A^(P, n, I, m) 



M^/n(n + 2) \ io n a 



n(n + 2) 



, <9 X , de, do 



(125) 



One may readily verify that the two transverse modes and the scalar Proca mode all have the property SA(X, n, I, m) = 
0, for A = M,E,P and therefore solve the Proca equation l|19fl . Moreover, these modes obey the (pseudo)- 
orthonormality conditions of Sect. IIVXI 

Since the Einstein spacetime is ultra-static, we can find the quantum weak energy inequality from Eq. (|109f) . To 
begin, let the worldline for a stationary observer be given by 7(t) = (r, x, 0, f) where (x, (9, ip) are the coordinates of a 
fixed point in space. An orthonormal basis in the neighbourhood around 7 is given by = \/|g' tM |(<9 / dx 1 ) 11 [no sum 
on /j]. All of the work is in the evaluation the classical energy density per mode, T(A, n, I, to), defined by Eq. (|107Jl . 
which, in the Einstein spacetime takes the form 



£(A,j) 



1 



MA,j) 



l-Foi I 



112 



o:>, 



1^12 



\T- 



13 



9-2 1 9-2 -2 

a z sin x a sm X sm 1 



a 4 sin 2 x a 4 sin 2 x sin 2 1 



a ■ 4 '2 

a 4 sm x sm 1 



-M 2 \Ao\ 



\A 2 



\A 3 



a 2 sin 2 x a 2 sin 2 x sin 2 1 



(126) 



Beginning with the ± modes, it is straightforward to calculate the field strength jF(±,n, l,m) = dA(±,n,l,m). 
Inserting the field strength into the above expression, and using d x ^ n i m — d x ^! n i m , we arrive at 



T(±, n, I, to) 



4a 4 r 



sin 2 x 



*(/ + !) 
sin 2 x 



1 



sin 2 x 



1 



sm 



\d X d^nlm\' 



. 2 - {d^nlml' 

sm 



(127) 



Note that the + and — modes have equal energy densities. Next, we use the definition Eq. (jll9[l of the scalar functions 
^nim and two identities. The first identity, 



1 1)1 I 



-At \d v Y lm \ 2 = -\A S 2 \Y lm \ 2 + 1(1 + 1) \Y lrn \ 2 
sm I 



(128) 
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is for the spherical harmonics where 



A.C2 = 



-da sin 9 da 



is the two-sphere Laplacian. The second identity is 



1(1 + 1) 

■ 2 

sin x 



Ki + (d x v nl y = u + -,K y*n 



—8 2 



1 



which follows from the ordinary differential equation satisfied by V n i (x) . Using both of these we arrive at 



X(±, n, I, m) 



2ah 



1(1 + 1) 
n l sin 2 x 



2n ' + -d' x )V^\Y lm \U 



1 



2 sin 2 x 



(d x V ni y + n z V^ 



^■S 2 \Ylm 



(129) 



(130) 



(131) 



Because the frequency uj n is independent of I and m, we can sum over these labels, obtaining considerable simpli- 
fication from the summation theorems for the spherical harmonics and Gegenbauer polynomials. In the case of the 
spherical harmonics it is well-known that (6l| 



E \ Y lm(0,V)t 



21 + 1 

47T 



(132) 



which reduces the sum over all / and m of T(=b, m) to 

7i— 1 +/ n — 1 



E E n±,n,is 



^7 E ^ I 2n 



1=1 m=—l 



1=1 



sin 2 x 



1 



»7ra J n sin x 

Using the summation theorem for the Gegenbauer polynomials 

n-l 



2~* 

n-l 



saw. 



nl • 



l-(2n 2 + k)sin 2 xE 



(21 + 1)1(1 + 1) 



i=i 



sin 2 x 



V, 



nl • 



E (2 ' + "'" +1 V .s M 



i=i 



[Eq. 8.934.3 of j£2|] we finally obtain 



sin 2 x 



n-l +2 



1 / sin n\ 



n \ smx 



E E 1(±,n,l,™>) = 



l—l m— — l 



(n 2 - l)u r , 
An 2 a? 



(133) 



(134) 



(135) 



It is noteworthy that this expression is exactly the total zero-point energy of n 2 — 1 harmonic oscillators at frequency 
Lo n , divided by the spatial volume 2n 2 a 3 of the Einstein universe. 

The contribution from the scalar Proca modes is calculated in much the same manner. The only difference is that 
we make use of the following two identities, 



(d x Tlm) 



1 fciln 2 , = „(„ + 2)n ni + ^ - - 

sin x 2 sin x 



A sin 2 x<9 x II 2 ; , 



(136) 



which follows from the ordinary differential equation satisfied by II n z and 

2(n + l) 2 



E(^ + i)n 2 i (x)- 



(137) 



1=0 

which again follows from the summation theorem for the Gegenbauer polynomials. For the Proca mode we find 

n +1 



E E Z(P,n,l,rn) = 



(n + 1) V, 
An 2 a 3 



(138) 
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Again we have found the zero-point energy times the multiplicity divided by the spatial volume. 

Finally, let s be the total number of spin degrees of freedom for the field, i.e. s = 2 for Maxwell and s = 3 for 
Proca. Substituting Eqs. Ijl35(l and Ijl38|l into Eq. IjlOBI) . we can write the quantum weak energy inequality in the 
Einstein spacetime for an arbitrary real- valued test function g(r) G C^°(R) as 



dr(:p:Ur)g(r) 2 > 



2 

n; | j 

(139) 

which holds for all Hadamard states u>. 

By suitably restricting the class of allowed states (if necessary) this bound may be extended to noncompactly 
supported g whose Fourier transforms have sufficiently rapid decay |75| . For electromagnetism, the above bound was 
evaluated for a Lorentzian test function in |27| ; here we will treat the example of a Gaussian test function in both the 
electromagnetic and Proca cases. Setting 

sM^-^r-^exp^^) , (140) 
the weight function <?(t) 2 is a normalised Gaussian. Using the Fourier transform of g, 

g(a) = 7r 1 /4(2r ) 1 / 2 ex P (- 



« 2 r 2 



(141) 



it is now possible to calculate 



da \g(a + u>)\ = 7rerfc(o;ro) 



(142) 



where erfc(cc) is the complementary error function. Accordingly, the quantum weak energy inequality for this test 
function is 



3s 

dr(: p :) w (t) 5 (t) 2 > - -^-j S s {t / a) 



(143) 



where the scale function 



- V (n 2 - l)\/n 2 + M 2 a 2 erfc( V / n 2 + M 2 a?z) 

Q ^ 



s- 2 



71=2 

oo 



J2( n + !) 2 y/n{n + 2) + M 2 a 2 edc(y/n{n + 2) + M 2 a 2 z) 



(144) 



is plotted in Figure For very short sampling times the scale function is close to unity and the QWEI bound is 
the same as that for Minkowski spacetime [22] . This makes sense because the spacetime appears to be flat over time 
scales which are short relative to the radius of the universe, and we should expect to recover the Minkowski space 
result. As the sampling time becomes progressively longer, the field has more time to "sample" the curvature of the 
universe, and thus it becomes increasingly difficult to generate negative energy densities. Also, as the mass of the 
field increases it becomes increasingly difficult to generate negative energy densities. This is seen by the faster decay 
in the scale functions for larger value of the field mass as seen in previous evaluations of the scale function for scalar 
fields in other spacetimes |15l Il7l Il8| . As mentioned above, the quantum inequality bound has also been evaluated 
for electromagnetism with a Lorentzian sampling function 27] and the behaviour is very much the same. 



APPENDIX A: CONSISTENCY OF THE HADAMARD CONDITION WITH THE CONSTRAINTS 



In this appendix we prove the following local-to-global result. 

Theorem A.l. Let u be a state on either 21 a/ (M,, g) or 2l(Al, g) and suppose that there exists a causal normal neigh- 
bourhood M of a Cauchy surface E in M. , and a Hadamard bisolution Wm ( respectively, W ) to □ + M 2 ( respectively, 
□ ) onAfxAf such that 



u (A(h)A(f 2 )) = W M (/i, (1 - M~ 2 d8) f 2 ) 



(Al) 
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for all fi S fig (AT) th e Proca case) or 

u([A](f 1 )[A}(h)) = W(f 1 J 2 ) (A2) 

for all fi E Oq(A/") ™£/i £/i = (in the Maxwell case). Then u> is Hadamard. 

Thus the Hadamard conditions introduced in Sect. IIII CI respect the Cauchy evolution of the field equations: if 
they are satisfied near one Cauchy surface, they hold near all others and hence hold globally. As the Hadamard 
form for Klein-Gordon bisolutions also propagates in this sense, our discussion may be regarded as checking the 
consistency of the definitions of Sect. IIII Cl with the Proca and Maxwell constraints. The propagation property is one 
of the key ingredients used in Sect. HVEl where we established the existence of Hadamard states on general globally 
hyperbolic spacetimes obeying our usual topological conditions. Although we will use the fact that S is compact for 
both the Maxwell and Proca fields this condition can probably be dropped; however in the Maxwell case it appears 
to be necessary [in Lemma |A.4f b)] to assume that -ffi(E) is trivial and that the compact support cohomology group 
H^(Ai) of the spacetime is therefore also trivial. 

The proof of Theorem I A . 1 1 relies on a number of results concerning the classical Proca, Maxwell and Klein-Gordon 
fields. One fact which will be used repeatedly is that, for / S SIq(M), E M f — if and only if / E (□ + A/ 2 )Oj(A4). 
Indeed, / = (□ + M 2 )E M f and since, by hypothesis, E^f = E M f, the support of E M f is contained in the compact 
set J + (supp/) n J~(supp/). Our first observation generalises this fact to the Proca field. 

Lemma A.2. Let M > 0. For f E Ql(M), A M f = if and only if f € (-<5d + M 2 )fiJ(A4). 

Proof. Since (1 — M~ 2 d8)(— 8d + M 2 ) = □ + M 2 , sufficiency holds because Em is a Klein-Gordon bisolution. 
Conversely, Am/ = is equivalent to Em(^ — M~ 2 dd)f = 0; hence, we have 

(-dS + M 2 )f = M 2 (a + M 2 )g (A3) 

for some g E f2g(A4). By applying Sd to both sides, we obtain Sdf = (□ + M 2 )6dg, which may be subtracted from 
the previous expression to yield 

(□ + M 2 )f = (□ + M 2 ){-5d + M 2 )g . (A4) 
Since both / and g are compactly supported, we conclude that / = (— Sd + M 2 )g. □ 

Proposition A. 3. (a) Given any f E Oj(A^) there exists g E £Iq(A4) such that f = f + (—Sd + M 2 )g is an element 

(b) Given any co-closed f E fl^Ai) there exists g E Q,q(A4) such that f = / — Sdg is a co-closed element ofCl^Af). 

Proof, (a) Choose smooth functions on Af with \ + + X~ = 1 an d X + equal to unity to the future of Af and 
vanishing to the past of J\f. Then 

f=(-8d + M 2 ) x + A M f (A5) 

belongs to Qq(AT) and satisfies A^/ = A^/. Applying Lemma rA.2l to / — /, it follows that / = / + (— 3d + M 2 )g 
as required. 

(b) With as above, set / = — Sdx + Eof. Then / belongs to f2g(AT) and is co-closed; moreover Eof and Eof are 
gauge equivalent by the proof of Prop. 4(c) in [3(|. Accordingly, there exists r\ E f2°(A / f) such that Eo(f — /) = dn. 
Taking co-derivatives, using SEq = EqS and co-closure of / and /, we have Sdr/ — and hence r\ — i?oC f° r some 
C E 0,q(M). Substituting back, we see that E (f — f — d£) = 0, so 

7-/-dC = n<? (A6) 

for some g E Qq(A4). Taking co-derivatives again, □(£ — dg) = and hence (because both £ and g are compactly 
supported) ( = Sg. Substituting back in Eq. (|A6|I . we have / = / + ddg + Ug = f — Sdg as required. □ 

Lemma A. 4. (a) If A is a weak one-form (□ + M 2 ) -solution obeying 

A{{-8d + M 2 )f) = {) (A7) 
for all f E HJ(A0, then Eq. holds for all / G Oj(A4). 

(b) If A is a weak one-form U-solution vanishing on co-closed f E CIq(AP), then A vanishes on all co-closed f E Qq(A4). 
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Proof, (a) Since, for / G f2j(A0, -4((n + M 2 )f) = 0, we have by Eq. JlTj that 

>4(d5/) = (A8) 

for all / G Qq(A/). Thus, .A(d-) is a global weak scalar (□ + Af 2 )-solution vanishing on 8£Iq(AT). Now we may fix 
fo £ fig (.AT) such that any / £ fig (A/ - ) may be written 

f = fojfdvol g + 6h (A9) 

for some h £ Oq(A/"). (Since Af is connected, boundaryless and orientable, this follows by de Rahm's theorem: see the 
remarks following Theorem 7.5.19 in Combining Eqs. (|A8I and (|A9|) . we have A(df) = „4(d/g) j f dvo\ g for all 

/ £ flo(Af). Accordingly, -4(d-) is constant on Af and hence (since it is a (□ + Af 2 )-solution) on A4. It follows that 
A(ddf) = A{df ) J Sf dvol g = for all / £ Qj(M). Since A((a + Af 2 )/) = 0, we deduce that Eq. Ezj) holds for all 
/ £ SIq(M) as required. 

(b) A(S-) is a two- form weak D-solution vanishing on all h £ fig (Af) and hence on all h £ Hq(A4). But since H^(A4) 
is trivial, any / £ Cl^(A4) such that 6f — can be written as / = Sh for some h £ Q,q(A4) and we conclude that 
A(f) = A(Sh) = 0. □ 

After these preliminaries, we may now prove the main result of this section. 
Proof of Theorem \A.1[ The arguments for the two theories run along parallel lines. First one extends the Wm 
(resp., W) to be a global bisolution to the appropriate one-form Klein-Gordon equation. Since the Hadamard form 
propagates, it suffices to show that Eq. IjAlfl holds for all /, £ ^(M) (resp., that Eq. (|A2|) holds for all co-closed 
/, £ £l\(A4)). In the Proca case, we apply Proposition lA.3f a) together with the field equation axiom P3 and Eq. l|AT|l 
to show that 

w (A(h)A{f 2 )) = lo (A(fx)A{h)) = W M (fx, (1 - M- 2 d6) f 2 ) (A10) 

for general ft £ Cll(A4). In the Maxwell case, we apply Proposition I A . 3f b) . axiom M3 and Eq. (|A2|) to obtain 

w([^](/i)M(/2)) =w ([AjifxMih)) = W(hJ 2 ) (All) 

for co-closed fi £ Hq(A4). The following claim will be proved below. 
Proposition A. 5. In the Proca case, 

W M ((Sd + M 2 )f, f) = W M (f, (Sd + Af 2 )/) = (A12) 

for all f,f £ Qq(A4), while in the Maxwell case, we have 

W(8df,f') = W(f',8df) = (A13) 

for all f, f £ fll(A4) with f co-closed. 

In combination with the explicit form of fi — fi (and, in the Proca case, the fact that (1 — Af ~ 2 d<5) and — 8d + Af 2 
commute) Prop. lA~5l allows us to show that 

W M (fx, (1 - Af- 2 d<5)/ 2 ) = W M (fi, (1 - Af- 2 d<5)/ 2 ) (A14) 

in the Proca case, and that W{f\, f2) — W{f\,fi) holds for the Maxwell field. Taken together with Eqs. (|A10|I 
and (|A1 l|l respectively, these relations then establish Eqs. IjAlfl and (|A2|) . □ 

The claim made above is proved as follows. 
Proof of Proposition \A.fA Fix an arbitrary /' £ J7q (Af) (with the additional requirement that Sf — in the 
Maxwell case). Then Wm(-,/') and !%(/', •) (or W(-,f) and W(f, ■) in the Maxwell case) obey the hypotheses 
of Lemma fA.4r a) owing to Eq. (|A1|) and axiom P3 (resp., Eq. (|A2|) and axiom M3). Accordingly, Eq. (|A12|I (resp., 
Eq. (|A"T3|n holds for all / £ Qj(M) and the fixed /' € ^o(^0' Now fix / G fll(M) arbitrarily. In the Proca case, 
Wm({— Sd + Af 2 )/, •) and Wm{-, (—Sd + Af 2 )/) are weak (□ + Af 2 )-solutions vanishing on Af and hence globally, 
as required. In the Maxwell case, W(Sdf, •) and W(-,Sdf) are weak O-solutions vanishing on co-closed elements 
of fl^Af). Lemma lA^f b) entails that they therefore vanish on all co-closed /' G ft})(A4), thereby completing the 
proof. □ 
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APPENDIX B: CONSTRUCTION OF A HADAMARD (□ + Af 2 )-BISOLUTION IN ULTRASTATIC 

SPACETIMES 



In this appendix we prove Theorem IIV. II bv explicitly constructing a Hadamard (□ + M 2 )-bisolution Wm on any 
ultrastatic spacetime (M.,g) obeying the assumptions stated at the beginning of Sect. II VI For M > 0, the argument 
proceeds as follows. First, we use functional calculus on the Hilbert space H to define Wm as a bilinear map from 
Hq(A4) x J7q(vW) to C. We show that Wm is in fact a one-form bidistributional weak (□ + M 2 )-bisolution, with 
antisymmetric part —iEm, and determine a crude bound on its wave-front set. Next, we appeal to the existence 
of a Hadamard (□ + M 2 )-bisolution Wm, also with antisymmetric part —iEm, established in Lemma 5.4(a) of [32| . 
A simple microlocal argument is used to show that Wm — Wm is smooth, from which it follows that Wm is itself 
Hadamard. Finally, we show that Wm may be expanded in terms of any K-pseudo-orthonormal complete set of 
eigenvectors for the operator K as claimed in Theorem IIV. II The argument is only slightly different in the case 
M = 0. 

It will be convenient to regard each A G fi 1 (A / t) as a smooth one-parameter family t A(t) of elements in 
H = A°(S) © A 1 (E), where A(t) is the restriction of A to the constant time surface {t} x S. The pairing (•, •) of 
one-forms on M. is related to the inner products of H and K by 



(A,B) = J dt((A(t) | B(f))) = Jdt (A(t),JB(t)j H 



(Bl) 



where J = 1a°(e) © -Ia^e)- We note that J and K = (-A| + M 2 ) © (-A s + M 2 ) commute. 



With these conventions, the Klein-Gordon equation (□ 
differential equation 



M )A = may be written as the Hilbert space ordinary 



~d 2 A(t) = KA{t) . 



(B2) 



Suppose that M > 0. Because E is compact, K has discrete spectrum bounded below by M 2 > 0. Accordingly, the 
operator K~ 1 / 2 is well-defined and bounded, and the advanced-minus- retarded solution operator Em '■ £Iq(M.) — > 
^(M) may be written 



{E M J){t) = / dt' A- 1 / 2 sin \K^ 2 (t' - tj\ Jit') , 



(B3) 



thus obtaining 



E M (J,J') = (J, E M J'} = / dtdt' (j(t), JK^ 2 sin [a 1 / 2 ^' - t)] ^"(0) H 



(B4) 



We define our candidate Hadamard (□ + M 2 )-bisolution Wm by taking the positive frequency part of —iEM, i.e., 
by replacing the sine function by an exponential. To be precise, for /, /' e D(M) and g, g' € fi°(£) © H 1 (E), we define 



W M {f®gJ'®g') = -\ QiK^K-^f'iK^Jg')^ , (B5) 

in which operators such as /'(A 1 / 2 ) are defined by functional calculus. Using the Cauchy-Schwarz inequality and 
elementary operator norm estimates, we find 

|WM(/®3,/'®<? , )l<^ll^- 1/2 ||||.g||||.g'|| sup |/(- w )| sup , (B6) 

from which it follows that Wm extends to a distribution in x S}q(A4))'. It is straightforward to check that 

Wm is a weak (□ + A/ 2 )-bisolution, with antisymmetric part —iE M - 

The wave-front set of Wm may be estimated in two ways. First, because it is a (□ + M 2 )-bisolution, we have 

WF {Wm) CJVxJV (B7) 

where M — {(x, k) e T* M. : g ab (x)k a ki, — 0} is the null bundle of (A4,g). Secondly, the explicit bound l|B6(l . coupled 
with the observation that sup wgo . ( - K i/2) |/(A — u>)\ is bounded as A — > +oo, but rapidly decaying for A — ► -co (and 
vice versa for sup IJEff ( K i/2) |/(A + u>)\) entails that 



WF (W M ) C (T+ U Z) x (T~ U Z) 



(B8) 
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where T ± = {(x, k) e T*M : ±k > 0} and Z = {(x, k) e T*M : k = 0} is the zero section of T*M. Comparing 
these two estimates, 

WF (W M ) C (7V+ U Z) x (Af- U Z) (B9) 

where 7V ± = TVn T ± are the future (+) and past (— ) null bundles. 

We now appeal to the existence of a Hadamard form (□ + M 2 )-bisolution Wm on (Ai, g) with antisymmetric part 
—iEm (Lemma 5.4(a) of Since the wave-front set of Wm also obeys (|B9|) . we have 

WF (W M - W M ) C (JV+ U Z) x (JV" U Z) . (BIO) 

But Wm — Wju is symmetric, so we also have 

WF (W M ~ W M ) C U Z) x (7V+ U Z) . (Bll) 

Comparing these two bounds, we see that WF (Wm — Wm) Q Z x Z and, since the wave-front set excludes the zero 
section, we conclude that this wave-front set is in fact empty. Accordingly, Wm — Wm (mod C°°) so Wm is of 
Hadamard form. 

The analysis of the massless case is complicated by the existence of a zero eigenvalue mode ipo — (1,0) for K. 
(Triviality of i?i(S) precludes the existence of any harmonic one-forms on S, so ipo is the unique zero mode.) However, 
the spectrum of K is otherwise bounded away from zero, so K~ x / 2 is well-defined and bounded on PH, where P is 
the orthogonal projector onto the orthogonal complement of <pq. In this case, we define 

W(f ®g,f®g') = -\ Q(K l '*)Pg, K-^fiK^PJg') ^ . (B12) 
It is easy to check that IF is a bidistributional D-bisolution, whose antisymmetric part is 



W(J, J 1 ) - W(J', J) = -iE (J, J') + i J dt dt' (t 1 - t) (J(t), ipoh (fo, JJ'(t')h (B13) 

and is therefore equal (mod C°°) to —iE . Appealing as before to the existence of a Hadamard D-bisolution W with 
antisymmetric part —iEo, the argument used above shows that W = W (mod C°°) because W — W is symmetric 
(mod C°°). It remains to show that the last term in Eq. 1|B13J| vanishes if J and J' are both co-closed, as required 
for consistency with the commutator axiom M4. Using Eq. (|B1(1 and the fact that Jipo = ipo, one may show that the 
term in question is proportional to (J, At) (tdt : J') — (J,tdt) (At, J'), which vanishes because (J~,dt) — (8J,t) = 
and similarly (dt, J') = 0. 

Finally, let £j be a complete set of K-pseudo-orthonormal iif-eigenfunctions, with corresponding eigenvalues io\ 
(u)\ > 0). Using the completeness relation Eq. I)43|l . we see that 

fi(KV*)*K-V*j(KV*)Pg = ]T ^f^J^)^ \ fc» \ g)) (B14) 

and hence 

W(J,J>) = - J2 ^m\tj))((J("i)9\m&\f'("j)9')) 

= - E ^r<fe i ( A i> J ) > (bis) 

where we have defined the modes Aj by Aj(t) = e~ luJjt !;\ and used 

| f>(u>i)g')) = J dte^dtj | f(t)g')) = J dt{{A 3 (t) \ J'(t))) = (A],j') (B16) 

and ((%j)g \ fc» = {J,A 3 ) = (Aj,J). 

In the case M > 0, exactly analogous arguments show that the Wm(J , J') may also be written in the form of the 
right-hand side of Eq. I|B15|I [the restriction to uij > is inessential as all modes obey this condition] . This completes 
the proof of Theorem lIV.ll 
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We also observe that their work is directed towards a construction of Hadamard states in general globally hyperbolic 
spacetimes which avoids the use of deformation arguments. 

Some care is needed over the definition of the two-point function in the Maxwell case (see Sec lllH . In particular, it is not 
a distribution. 

The support of a function is the closure of the set of points on which it is nonzero. 

The index is the number of spacelike (i.e., negative norm-squared) basis vectors in any gr-orthonormal frame. 
The topology is defined so that a sequence U n — > in fiJ^A/") if and only if *(V A *U n ) — > in T>(J\f) for every V € Q, p (Af). 
Since £ is compact, it admits a finite good cover [Theorem 5.1 in |63fl ] which determines a finite good cover of M = R x S. 
Accordingly, we may use Poincare duality [Theorem 5.4 in |63ll ] to prove that H^(M) is isomorphic to H 1 (M) and hence 
to i? 1 (E) [see the remarks following Corollary 5.1.1 in [(jUll. This is isomorphic to .ffi(E) by de Rahm's theorem [64[. 
Note that Dimock |36l| uses E for the retarded-minus-advanced bisolution. 

For example, 21m (A^, g) is the infinitesimal Weyl algebra over the real symplectic space Qo(A4; R)/ ker Am (with the 
quotient topology) with the addition of the *-operation. This has no nontrivial two-sided *-ideals (see Sect. 7.1 in |4(J). 
Effectively, one is substituting the test function f(x) = J dr\h{x)\~ 1 ^ 2 5 A {x — 7(i~))<7 Q (r) in Eq. (1761 . This results in the 
required expression. 

The RHS of the bound may be written in the form — J °° duQ(u)[g(u)\ 2 , where Q(u) is a polynomially bounded function. 
Accordingly, the integral will converge provided that the Fourier transfor m g (u) decays sufficiently rapidly as u — > +oo; 
equivalently, provided that g belongs to a Sobolev space VK fc,2 (R) (see e.g., |6fl| 1 of sufficiently high regularity (k sufficiently 
large). If we restrict the class of states so that (: p :)^(t) belongs to the dual Sobolev space W- k ' 2 (M) then Eq. 
will continue to hold. To see this, we approximate g in the W k ' 2 (M.) norm by a sequence of compactly supported smooth 
functions [using density of Co°(R) in W ,2 (R)] and apply Eq. 113911 . It is not presently known whether this requires a 
nontrivial restriction on the class of allowed states. 
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FIG. 1: Plot of the scale function S s (z) for the Gaussian sampling function in the Einstein Universe. Three different cases are 
of Eq. (11441 are displayed: Electromagnetism (solid line) and two cases of the Proca field for different values of the mass times 
the radius of universe, Ma = 1 (dashed line) and Ma — 2 (dotted line). The plots were produced by evaluating the first 3000 
terms in the summation of Eq. 11441 . 



